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Using Lie group theory and canonical transformations, we construct explicit solutions of nonlinear 
Schrodinger equations with spatially inhomogeneous nonlinearities. We present the general theory, 
QQ I use it to study different examples and use the qualitative theory of dynamical systems to obtain 

, some properties of these solutions. 

o 

■ PACS numbers: 35Q51, 35Q55, 34C14. 
0\ \ I. INTRODUCTION 

^0 . The Nonlinear Schrodinger Equation (NLSE) in its many versions is one of the most important models of mathe- 
Ph ' matical physics, with applications to different fields ^33*1 as for example in semiconductor electronics [1, [2l|, nonlinear 
, optics [IS] , photonics [13 , plasma physics [fO] , fundamentation of quantum mechanics [1^ , dynamics of accelerators 
[ist . mean- field theory of Bose-Einstein condensates [1, [HI or biomolecule dynamics 0] to cite only a few examples. 
In some of these fields and many others, the NLS E ap pears as an asymptotic limit for a slowly varying dispersive 
wave envelope propagating in a nonlinear medium [30|. Moreover, the range of applicability is large because of the 
well-known universality of this equation 
^ , The study of these equations has served as a catalyzer of the development of new ideas or even mathematical 
J — ' concepts such as solitons [S^ or singularities in partial differential equations [l^, |3l| . 

In the last years there has been an increased interest in a one-dimensional nonlinear Schrodinger equation with 
inhomogeneous nonlinearity (INLSE): 

I ■ iipt = + V{x)^ + g{x)\il}\^ i^, (1) 

o ; 

00 , with a; G R, V{x) is an external potential and g{x) describes the spatial modulation of the nonlinearity. This equation 

■ arises in different physical contexts such as nonlinear optics and dynamics of Bose-Einstein condensates. But it is in 
' the later field where the possibility of usiiig the Feschbach resonance management techniques to modify spatially the 

• ^ . collisional interactions between atoms [11, [Til, [IB, [IE [H, [13, [13, [3 the one which has motivated a lot of theoretical 

■ research in the last few years focusing on questions of direct applicability to experiments. Different aspects of the 
^ , dynamics of solitons in these contexts have been studied such as the emission of solitons |27|.|33 | and the propagation 

of solitons when the space modulation of the nonlinearity is a random T|, periodic (25l. |29||. linear [s^l or localized 
function [26]. Although many exact solutions of the NLSE with spatially homogeneous nonlinearities and without 
potentials {V = 0) have been known for a long time, the problem of finding exact solutions even of the NLS with 
homogeneous nonlinearities and general potentials is a very difficult one. 

In this paper, using Lie symmetries we find general classes of potentials V{x) and nonlinearity functions g{x) for 
which exact solutions can be constructed by combining solutions of the integrable NLS and solvable potentials V{x). 
The basic idea of the Lie symmetries method is to study the invariance properties of given differential equations under 
continuous groups of transformations. This method has been applied successfully to different equations, such as, for 
example, differential equations that model anharmonic oscillators [l^ 113] and Madelung fluid equations [s^] . In Ref. 
[3] we have presented some examples of this methodology of specific physical interest. Here we complement that 
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analysis by presenting the general theory, provide more examples, study the case of asymmetric solutions and use 
qualitative theory of dynamical systems to provide a much more complete analysis of the method and its applications 
to equations of physical relevance. 

The paper is organized as follows. In Section 2, wc introduce the general theory of the Lie symmetry analysis 
for ordinary differential equations and particularize it for our model problem: the nonlinear Schrodinger equation 
with an inhomogeneous nonlinearity. In Section 3, we study the canonical transformations of the INLSE. In Section 
4, we present the connection between the NLSE and the INLSE. In Section 5, we use the method to construct 
explicit solutions of the stationary nonlinear Schrodinger equation with an inhomogeneous nonlinearity and study the 
qualitative behaviour of the NLSE and its qualitative connection with the INLSE in different examples. Finally, in 
Section 6, we present asymmetric solutions of the INLSE. To our knowledge, this is the first time that such solutions 
are calculated. 



II. GENERAL THEORY OF LIE SYMMETRIES 



In this paper we will look for localized stationary solutions of Eq. ([T]), which arc of the form 

V'(a;,t) = u{x)e-'^\ 

which satisfy the following nonlinear eigenvalue problem 

^ + V{x)u + g{x)u^ — Au, 
lim u(x) — 0. 

|a:| — >oo 

By definition 0, H^Ij ^ second-order differential equation A{x,u,u' ,u") — possesses a Lie group of point transfor- 
mations or Lie point symmetry of the form 

M = £_{x, u)d/dx -t- ?7(x, u)d/du, 

if the action of the second extension of M, i.e. M^^^ on A is equal to zero, i.e. 



(2a) 
(2b) 



M(2)A(a;,'u,u',u") 



with Ty'*^' given by 



(■\ \ , \ ^ (2) / \ ^ 



A{x, u, u', u") = 0, 



{x,u,u',u", ...,1*") = — u'^'^' ^ \ 



Dx 



Dx 



where ry*-*^^ — ri{x,u) and D/Dx is the total derivative, i.e. 



D d , d 



d 



Dx dx di 



U — h U 



For example, rj'^^^ is equal to 



In our case, A(x,u,Ux,Uxx) is given by 



A{x, U, Ux, Uxx) — ~Uxx + fi^i u), 



where f{x, u) — V{x)u + g{x)u^ — Am, and the action of the operator M^^^ on A{x, u, Ux, Uxx) leads to a polynomial 
equation in Ux- By equating coefficients of powers of Ux, one obtains 



'^Vxu ^xx '^f^u — 0: 
Vxx - ?/x - Vfu + Vuf - 2^xf = 0. 



(3a) 
(3b) 
(3c) 
(3d) 
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Integrating Eqs. ((3a|) and pb| . we get 

S,{x , u) — a{x)u + b{x) , ri{x , u) ^ a' {x)u'^ (x) + c{x)u + d{x) . (4) 
Substituting these expressions into Eq. ([5c|) we obtain 

2c'{x) = b"{x), a(x) = 0. (5) 
Finally, substituting Eqs. ^ and (O in Eq. (pd]). we get 

^(XjM) = &(a;), (6a) 

ri[x,u) = c{x)u, (6b) 

c"{x)~b{x)V'{x)~2b'{x){V{x)- X) = 0, (6c) 

2c(a:)5(x) + fe(.T),g'(x) + 26'(.T)g(x) = 0. (6d) 

Also the substitution of Eq. (O in Eq. (j6dp gives 

g(x) =go&-'(x)e-'^^o"=i/K^)<i«. 

where go and C are arbitrary constants. Summarizing the previous calculations, the Lie point symmetry is of the 
form 

M = 6(.)|- + c(x)u|-, (7) 

where 

g{x) ^ gob-\x)e-'^fo'Ms)ds^ (8^) 

c{x) = ^b'{x) + C, (8b) 
c"{x)-b{x)V'(x)-2b'{x){V{x)-X) = 0. (8c) 

Eqs. ([8]) allow us to construct pairs {V{x), g{x)} for which a Lie point symmetry exists. Thus given either g{x) or 
V{x), in principle we can choose the other in order to satisfy Eqs. Q. In what follows, we will study the implications 
of the existence of this Lie symmetry. 



III. CANONICAL TRANSFORMATIONS AND INVARIANTS. 

It is known [loj . that the invariance of the energy is associated to the translational invariance. The generator of 
such a transformation is of the form M = d/dX. To use this fact, we define the transformation from variables (XjM) 
to new variables (X, U) 

X = h{x), U = n{x)u, (9) 

where h{x) and n{x) will be determined by requiring that a conservation law of energy type M = d/dX exists in the 
canonical variables. In fact, using Eq. ([9]), we get 

|:^„'(x)„^ + .'(x)^. (11) 

Inserting the expressions pUj) and (jlip in Eq. ^ and assuming the condition M = d/dX, one finds 

h'{x)b{x) = 1, 

b{x)n'{x) + c{x)n{x) = 0. (12) 
By inserting Eq. (I8bp into (fT2|) . and integrating, we obtain 

Jo Ks) 

n{x) = 6(x)-i/2g-c/„-i/Ks)ds_ 
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We can now write Eq. (|2a|) in terms of the canonical coordinates U and X, 



with 



E^{X~ Vix)) h{xf - \b\xf + U{x)h"{x) + C^. 



(13) 



(14) 



Equation (fT3|) is the so-called Duffing equation which arises as a model of damped nonlinear oscillations 14]. It follows 
from Eqs. (j8bp - ([5c|) that the quantity E given by Eq. is a constant of motion. 

When C — Q the previous transformations preserve the Hamiltonian structure, because the canonical transformation 
is symplectic. In that case, Eq. (fT5)) becomes 



dX2 



goU-^ = EU. 



(15) 



As E is constant, this means that in the new variables we obtain the nonlinear Schrddinger equation (NLSE) without 
external potential and with an homogeneous nonlinearity. 

Of course not all choices of V{x) and g{x) lead to the existence of a Lie symmetry or an appropriate canonical 
transformation, since they are linked by Eqs. ([8]). This fact imposes some obvious restrictions, for instance b{x) must 
be smooth and positive. 

Note that Eq. (fT5|) is a stationary homogeneous NLSE without external potential, which can be reduced to 
quadratures and for which many solutions are known. So, we obtain 



X-X^ 



dG 



'Uo y'2(iV + i^G2 + igoG4) 
with TV a constant of integration. Moreover, the energy of the system is given by 



H = 



1 f dU 



1 9 1 4 



2 \dX ^ 

Many solutions of Eq. (|15p are known. In this paper we will use the following ones 

1 



U,{X) 

U2{X) 

Ui{X) 
Ui{X) 



= V 



cosh{fiX) ' 

rj tanh(/iX), 

sn(/iX, k) 
^ dn{iiX, k) ' 

rj dn(/iX, k), 



E = -fi ,go = 2" 



9 2u2 
E = 2fi\go = ^ 



E^fi^l^2k^),go^ 
E = fi^ik^-2),go = 



2^i^k^{l-e 



2fi' 



(16) 

(17a) 
(17b) 
(17c) 
(17d) 



with < fc < 1. Table U summarizes the parameter values required for the existence of the solutions listed in Eqs. 



u 


E 


90 


H 


Ui{X) 


negative 


negative 





U2{X) 


positive 


positive 


positive 


UaiX) 


both 


negative 


positive 


Ui{X) 


negative 


negative 


negative 



TABLE I: Conditions on the parameters E, go and energy H for the existence of the solutions Ui,i = 1..4 of Eq. (|15|l listed 
in Eqs. ([T7)l . 
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IV. CONNECTION BETWEEN THE NLSE AND INLSE VIA THE LSE 

Setting C — and eliminating c{x) in Eqs. jS]) we get 

gix) = go/6(x)3, (18) 

plus an equation relating b(x) and V{x) 

b"'{x) - 2h{x)V'{x) + 46'(a;)A - Ab'{x)V{x) = 0. (19) 

We notice that the simplest form to generate solutions for our problem, which involves constructing solutions pairs 
{b{x),V{x)) of Eq. p9|) is to fix b{x) and then, to calculate V{x), since then we must solve a linear first order 
equation. Although we can eliminate b{x) and obtain a nonlinear equation for the pairs g{x) and V{x) for which there 
is a Lie symmetry, it is more convenient to work with (1 191). which is a linear equation. Alternatively, we can define 
p{x) = 6^/^(a;) and get an Ermakov-Pinney equation [13. |24| 

p,, + {X-V{x))p = E/p^ 

whose solutions can be constructed as 

p = {aipl + 2(3ipiip2 + IvlY'"^ , 
with a,P,"f constant and (pj{x) being two linearly independent solutions of the Schrodinger equation 

Vxx + (A - V{x)) ip ^0. 

This choice leads to £' = AW"^ with A ~ a'j — (3"^ and W being the (constant) Wronskian W = (p'i(p2 — ViV^- 
Thus, given any arbitrary solution of the linear Schrodinger equation (I20p we can construct solutions of the nonlinear 
spatially inhomogeneous problem Eq. Ii2a\) from the known solutions of Eq. Thus, using the huge amount of 

knowledge on the linear Schrodinger equation we can get potentials V{x) for which ipi and are known and construct 
6(x), the canonical transformations h{x),n{x), the nonlinearity g{x) and the explicit solutions u{x). 

V. QUALITATIVE ANALYSIS AND EXACT SOLUTIONS 

In this section, we will calculate exact solutions of Eq. ([^al) for different specific choices of the nonlinear coefficient 
g{x) and the external potential V{x), using the method described in the previous sections, for C — 0. Moreover, using 
qualitative analysis, we will describe properties of the solutions of Eq. (j2ap . on the basis of the qualitative behaviour 
of Eq. (fTS]) . We begin by calculating the equilibrium points of this equation. One easily finds that the Eq. (|15p has 
three possible equilibrium points, depending of the signs of E and go: 

U± = ±^/eI^o, 
U ^ 0. 

Then, we distinguish four cases: 

1. For E < 0, go < 0, we get three equilibrium points. U = is a saddle point and U± are centers. Fig. [Ua). 

2. When E < 0, go > 0, we obtain that U — is the only equilibrium point, which is a saddle point. Fig. [IKb). 

3. When E > 0, go > 0, we get three equilibrium points. U± are saddle points and J7 = is a center, Fig. [ijc). 

4. The last case corresponds to i? > 0, go < 0. For this case, the only equilibrium point is the trivial solution 
U — 0, which is a global center. Fig. [IJd). 

Using Eq. (fTO)) . we can draw the phase portrait of Eq. (fT5|) . as we can sec in Fig. [T] 
In what follows we will present three examples as applications of our theory: 

Example 1. Let us take b{x) = cosh(2;). By using Eqs. (fTO|) and p4p . for C = 0, we obtain 

=x + \+(\-e] 

^ ^ 4 V4 / cosh2(a;) 



dU/dX 




dU/dX 



dU/dX 




FIG. 1: [Color online] Phase portrait of the real solutions of Eq. (fTS)) for (a) E < 0, go < 0, (b) iJ < 0, go > (c) i? > 0, 
go > and (d) E > 0, go < 





FIG. 2: [Color online] Solutions of (a) Eq. ^ and (b) Eq. ((2al), for £ = 0.15, go ^ ~1 (solid blue line) and E = -0.75, 
go = — 1 (dashed red line) in both cases. We apply the transformation 1)21^ to the solutions of Eq. (|15|l shown in Fig. [IJa) to 
obtain the solutions of Eq. (|2ap shown in Fig [2f b) . 



Moreover, using Eq. (|18p . g(a;) is given by 



with X{x) being 



90 



cosh'^ (x) ' 



cosX{x) = — tanha;. 



where < X < tt, subject to the Dirichlet boundary conditions U{Q) = U{tt) = 0. Any solution U of Eq. gives 
a solution 



u{x) = b^^^{x)U{X{x)), 
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FIG. 3: [Color online] Example black solitons solutions of (a) Eq. (fTS)) and (b) Eq. (f2a)) with go = 1, X = 1/4 and (i) a = 0.1 
(dashed-dot green line), (ii) a = 0.4 (solid blue line) and (iii) a = 0.7 (dashed red line). The solutions shown in Fig. EJb) are 
obtained from those shown in Fig. IHa) through the transformations (|9]). 



of the original equation (|2ap . If £^ < and 170 < 0, we are in the first case, Fig. [Ija). The periodic solution (|17cp of 
Eq. (fT5|) is a closed orbit of the phase portrait ([/, dU/dX), corresponding to one of the external closed orbits to the 
homoclinic orbits, shown in Fig. [Ija). By an elementary application of L'Hopital rule, it is easy to verify that the 
solution 

u{x) = b'/^ix)U,iX{x)), (21) 



is a homoclinic orbit to zero (bright soliton) of the original equation ([2a|) . 

If -E > and go < 0, we are in the four case, where J7 = is a center, Fig. \V[d). Again, u{x) = b^^^{x)U3{X{x)) 
is a homoclinic orbit (bright soliton). The solutions of Eq. (|2ap for E = 0.15, go = —1 and E — —0.75, go — —1 are 
drawn in Fig. [2](b). The case E — 1/4, where V{x) is a constant, was studied in 

Example 2. Let us take V{x) = 0. Then Eq. becomes 6"'(x) + AXb'{x) = 0. For A > 0, the solution can be 
written as 

b{x) = 1 + acos{2\^x). 
Using Eq. psp . we obtain a periodic nonlinearity 

g{x) = go{l + acos 
For small a, this nonlinearity is approximately harmonic 

g{x) ~ go{l — iacos{2\^x)), a ^ 1. 
We can construct our canonical transformation by using Eqs. dH) and obt am 



X(x) = ^ arctan ( \ - — - tan(\/Aa;) | . (22) 

^ ^ v/A(l-a2) \ \ 1 + a ^ 'J ^ ' 

Using any solution of Eq. p5|) with E — X (^1 — a^^ this transformation provides solutions of Eq. ([2aP with g{x) given 
by ((22)) . For example, when go > we can use U2 as defined by (jl7bp . which in the phase portrait is the heteroclinic 
orbit shown in Fig. [TJc). So, the solution of Eq. ([2aP is of the form 



u{x) 



— ^ '- [l + a cos(2VAa;) x tanh \ — ^ ^ X(a;) 

.go V / V 2 



(23) 



As we can see in Fig. EJb), u{x) is a heteroclinic connection between periodic solutions of Eq. (|2ap. It is important 
to note that in the asymptotic regions the profile of u{x) is close to b^^^{x) multiplied by a constant. Therefore, the 
canonical transformation ([5]), in this case, transforms a heteroclinic orbit in the phase portrait (U,dU/dX) into a 



8 



300 





FIG. 4: [Color online] (a) Quasi-harmonic potential for M — 
red line), (b) Solutions of Eq. ([2a)l for a = 1, go = — 1, 
line) 



1, A = 1 and (i) /3 = 0.5 (solid blue line) and (ii) f3 = 2.5 (dashed 
= 1 and (i) (3 = 0.5 (solid blue line) and (ii) /3 = 2.5 (dashed red 



heteroclinic connection, Eq. ((23|l . On the other hand, any closed orbit U inside the heterochnic loop of the phase 
portrait (U,dU/dX) provides a new heterochnic connection of the original equation (|2ap . 

Example 3. The last example is the so-called quasi-harmonic confinement V{x) ~ x^. 
If we choose b{x) — a/ ^/l + f3x^, with a,/3 > 0, then, we obtain 

.90 



and 



V{x) = M{1 + I3x^) + 



4 {T+J^ 



with M a positive constant. Although the expression of V{x) is complicated this potential V{x) is a quasi-harmonic 
potential and satisfies V{x) ~ for large x. Moreover V{x) is a harmonic potential with a bounded perturbative 
term (see Fig. IDJa)). As to the nonlinear term, it satisfies, g{x) ~ x'^ for \x\ ^ 1, and g{x) ~ x^ for |a;| ^ 1. Using 
Eq. (|14p we get E — ^a^M . Taking < 0, we obtain the nonlinear Schrodinger equation with nonlinear attractive 
term, Eq. p5)) . As i? < and go < 0, all the solutions of Eq. p5)) are bounded, (see Fig. [Ija)). If U{X) is one of 
these solutions, it is clear that 

u{x) = b{xf'^U{X{x)), 

is a homoclinic orbit (bright soliton) of the original equation. In particular, the solution given by Eq. (|17aD is 



Ui{X) = 




1 



50 cosh{^/\E\X) 



As X{x) = xy^l + (3x^/{2a) + smh~\y^x)/{2ay^), we get 



uix)^b{xy/^UiiX{x)). 



In Fig. m^b), we draw the solutions of Eq. (Pa|l for different values of the parameter (3. 

In this example, we have used the solution Ui{X) of Eq. (fT5|) . Another possibility is to choose the closed periodic 
orbits inside the homoclinic loop (see Fig. [TJ^a)) given by {/^{X) in Eq. (|17dp . We note that, for this case, Eq. (fTB)) 
satisfies H <0, as one can see in Table HI 

The analytical expression of the closed periodic orbits outside the homoclinic loop is given by[/3(X), withfc > 1/V2. 
Such orbits satisfy H > (see Table |l|. 

As X is a bijective map on the real line, b{x) is a positive function and U3{X) is a periodic function with infinite 
nodes on the real line, the function u{x) = 6(a:)^/^J73(Ar) also has infinite nodes on the real line, as it is shown in Fig. 
[51 It is immediate to check that the solution u — > when x ±oo. Moreover, the zeroes of (j) accumulate for large 
values of x, since the distance between two consecutive zeroes is given by Xn+i — Xn ^ \/n + 1 — ^/n. Other localized 
solutions with infinite nodes have been studied in a different context in Ref . Q . 



-3 

-10 r. 
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FIG. 5: Solution of Eq. \2'a\ with infinite nodes for a = 1, /3 = 2.5, go = —1 and = 3/4 

VI. ASYMMETRIC MODES OF THE INLSE 

In the previous section, we have explored the case C = in the canonical transformation ([5]). In this case, the 
original equation ((2a|) is reduced to a Nonlinear Schrodinger Equation. In this section, we will study the case C > 0. 
If we take <Q and £' > 0, the resulting equation is 

^ + 2C- + MC/^ + i^f/ = 0. (24) 

In general, this equation is not integrable and the energy is not a conserved quantity. However, exact solutions of Eq. 
(|24p can be constructed analytically in particular cases. An exact integrability condition was given in [23 | 

In that case a family of exact analytical solutions of Eq. ([24| is given by the expression 

^"V2M dn(i^(l-e---),V2/2)' ""^'^^S,... 

where and B are related to the boundary conditions of the problem . 

We are going to solve Eq. ([2a|l. using the solutions of Eq. (p4|) . Choosing b{x) = cosh(a;) and using Eqs. ((T9|) and 
we can calculate the potential V{x): 

The nonlinear term is 

5(a:)=5oCosh-3(a;)e-2^^(-). 

So, we can calculate the solutions of Eq. ([2a|l for the case C ^ and to compare them with the solution obtained for 
the case C = 0, example 1. In this way, we can construct our canonical transformation by using Eqs. ^ and obtain 

cosA'(a;) = — tanhx. 

Then, < X < tt and using the boundary conditions for u, lim|j.|^oo u{x) — 0, one has to impose U{Q) — U{'k) — 0. 
Using these boundary conditions, we obtain the value of the amplitude jin as a function of the value of the integer 

n 

_ 4Cg(v^/2) 
- 3(l-e-2CV3)"' 



and B = 2C/3. Thus, the solutions of Eq. ^ are 



^-2cx/3f^ {2nK{V2/2)il - 6-2^)7(1 - e-^^-/^),V2/2) 



y2M dn {2nK{V2/2){l ~ e-2CX/3)/(l - e-2Cx/3), V2/2) ' 

n = 1,2,3,... 
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FIG. 6: Asymmetric solutions of Eq. (|2a|) with C — 4 for (a) n = 1, (b) n = 2 and (c) n = 3 



20 



where K{k) is the elHptic integral of the first kind, 

K{k)^ 



r/2 



de 



1 - fc2 sin^(6l) 



Then, the solutions of Eq. (l2a|) are 



^„(a;) =6i/2(^)gCX(x)^^(^(^))^ n=l,2,... (25) 

By using L'Hopital's rule, ^ when |a;| ^ oo in (|25p . So, the solutions (|25p are localized solutions of our 

problem as it can be seen in Fig. [6l These solutions are asymmetric solutions of Eq. ([2a|) . Moreover, each of those 
solutions has exactly n — 1 zeroes. In Fig. |6l we plot some of them corresponding to n = 1, 2, 3. The picture in Fig. 
[SJa) shows clearly the difference between the positive solution given by (pS)) . for C 7^ 0, and the positive solution 
plotted in Fig. ^h) and given by Eq. jH]), for C = 0. 



VII. CONCLUSIONS 



In this paper, we have used the method of Lie symmetries to find exact solutions of the INLSE. We have introduced 
the general framework of the Lie's theory and presented different examples as application to the theory. By using the 
qualitative theory of the dynamical systems, we can show the properties of the solutions of the INLSE and to classify 
such solutions. Finally, we have calculated asymetric solitons of the inhomogeneous nonlinear Schrodinger equation. 
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